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ABSTRACT 

It is predicted that the Ha emission hne at the early nebular epoch of type II-P 
supernovae may display robust observational effects of radiation damping wings. Monte 
Carlo simulations demonstrate that the Ha line with the optical depth typical of the 
early nebular epoch acquires extended wings and redshift and becomes broader com- 
pared to the Sobolev approximation case. The strength of these effects may be used for 
constraining parameters of the line-emitting zone. The anomalous redshift, width and 
red wing of the Ha emission in the supernova SN 1997D on day 150 are explained in 
terms of damping wings effects. 

Subject headings: line: profile, formation — supernovae: individual (SN 1997D) 



1. Introduction 

The ejecta structure of a type II-P supernova 
(SN II-P) may be sketchy portrayed as a combina- 
tion of a hydrogen envelope and a " mixed core" . The 
latter, being a macroscopic mixture of newly synthe- 
sised elements, such as. He, O, ^^Ni with a fraction of 
the hydrogen envelope, is responsible for the luminos- 
ity at the nebular epoch, i.e., at the age > 100 d. The 
model of mixed core of SN II-P is fundamental to the 
assessement of the nuclesynthesis yields, which is the 
subject of great interest. Two recent nebu lar models 
of SN 1987A (De Kool, Li & McCray [199^; Kozma & 



Fransson 1998) demonstrate that the basic physical 
processes in the mixed core are well understood to 
make feasible the composition analysis. However, to 
achieve this goal a nebular model of SN II-P should, 
first, pass the rigorous test on reproducing fluxes and 
profiles of hydrogen lines, including Ha. 

Here I argue that the model of the Ha line in 
SN II-P at the early nebular epoch (100-200 d af- 
ter the outburst) should essentially include effects of 
radiation damping wings always ignored before. The 
physical reason for the damping wings to play a role 
is a high degree of a non-thermal excitation of hy- 
drogen in the mixed core of SN II-P, which results 
in the large Ha optical depth at the early nebular 
epoch. In this situation to model Ha damping wing 
effects the Sobolev (local escape) approximation must 
be abandoned. The radiation transfer in the physical 
and the frequency space simulated with the Monte 
Carlo method is applied to demonstrate expected ef- 
fects of the Ha damping wings in SN II-P (Section 2). 
Recently found in SN 1997D a disparity between the 
large asymmetry of Ha at the early nebular phase and 
the symmetric Ha profile at the late nebular epoch 
(Chugai & Utrobin 1999) admits a simple explanation 
in terms of the damping wings (Section 3). Another 
motivation to address the issue of the Ha damping 
wings in SN II-P is that this factor may provide an 
independent diagnostic tool for constraining parame- 
ters of the mixed core. 

2. Expected effects of Ha damping wings in 
SN II-P 

The Sobolev approximation widely used to analyse 
supernova optical spectra is justified for most lines 
because of a highly supersonic expansion, v/u » 1, 
{v and u being the expansion and the thermal ve- 
locity, respectively). However, if the optical depth 



in damping wings is large, the "supersonic" criterion 
for the local escape of photons should be replaced by 
the inequality v/u 3> ar (Chugai 1980), where a is 
the Voigt parameter, r = k\t is the Sobolev optical 
depth {k is the line integrated absorption coefficient, 
A is the wavelength, t is the expansion time) . The Ha 
Voigt parameter is a = A21)/ {At:1^vd) = 

3.27 X 10"'^, where Aik is the Einstein transition prob- 
ability, Aj/d is the thermal Doppler width (kinetic 
temperature of 5000 K is adopted). At the early neb- 
ular epoch (~ 150 d) the Sobolev optical depth of the 
mixed core in Ha may exceed 10^ (Xu et al. 1992) 



implying that v/u ~ 2 x 10 , i.e., is of the order of 
or. In this case the local escape approximation is too 
crude, so that the full radiation transfer should be 
applied to compute the Ha line profile. 

One needs, first, to specify the frequency redis- 
tribution function for the Ha damping wings. Let 
x' and X be the incident and scattered photon fre- 
quencies in the atom frame. With lower and up- 
per levels broadened by the radiation damping, the 
atomic angle- averaged frequency redistribution func- 
tion in far wings {\x'\ ^ a) may be written as 
r(x' ,x) = 4>{x' )p{x' , x) (Oxenius 1986 ), where 4>{x) = 
{a/'K)/{x^ + a^) is the normalized absorption coeffi- 
cient and 
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Here 02 = A21 / {Att Iw^,) and 03 (A3i+y432)/(47rAt'D) 
are Voigt parameters for the lower and the upper 
level. Thermal velocity effects in the wings are neg- 
ligible and will be ignored here. This implies that 
the frequency redistribution function in the wings is 
the same for the atom and the fluid frame. The re- 
distribution function may be essentially simplified to 
reduce the number of scattering in the Monte Carlo 
computations preserving good precision. Since multi- 
ple local scattering in the Doppler core {\x\ < Xc ~ ^) 
eventually ends up with the photon re-emitted in the 
wings, the scatterings in the Doppler core may be 
omitted by adopting the function p{x' , x) equal to 
zero in the core |a;| < x^. After the re- normalization 
this function is labeled P{x',x). For instance, the 
function (j){x) in the second term of Eq. (|^) transforms 
into ip2{x) = 0.5xc/x^ for > Xc and ip2ix) = oth- 
erwise. Results are not sensitive to the choice of the 
core width and = 3 is assumed. 

The Ha photon history is simulated as follows. 
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The photon is launched with the distribution function 
il)2{x). When propagating it experiences a redshift in 
the co-moving frame until it randomly scatters or es- 
capes. If scattered, the frequency of the re-emitted 
photon is randomly choosen according to the modified 
distribution function P{x' ,x). The process continues 
until the photon escapes. 

To demonstrate the role of Ha damping wings in 
SN II-P let us consider a freely expanding envelope 
(v — r/t, V being the velocity at the radius r at the 
expansion time t) with a homogeneous line-emitting 
core {v < vi) and a scattering halo {vi < v < V2). 
The Sobolev optical depth of the halo decreases out- 
ward as r = T2(l — V /v2)- The core, optically thin in 
the continuum, may emit a continuum radiation from 
homogeneously distributed sources. The velocity and 
the optical depth in the core (subscript 1) and in the 
halo (subscript 2), respectively, are given in Table 1. 

The homogeneous core with the low optical depth 
in the model Ml produces a parabolic profile (Fig. 
la) expected for the Sobolev approximation. When 
the optical depth is high (models M2 and M3) the 
extended wings emerge, while the line core becomes 
broader and gets redshifted. All the effects are caused 
by the radiation transfer in the damping wings. Note, 
that in the case of the Sobolev approximation the 
models M2 and MS would have the same parabolic 
profile as that of the model Ml. A conspicuous kink 
at the blue slope in the model M3 is caused by a self- 
absorption due to non-local scattering. The models 
M4 and M5 show how the models Ml and M3 are 
changed when a scattering halo and a continuum core 
are present (Fig. lb). The model M5 gives a realistic 
impression of damping wings effects in Ha at the early 
nebular epoch of SN II-P. 

How would the inhomogeneity of the hydrogen dis- 
tribution in the mixed core change these results? 
Generally, the optical depth in the damping wing 
depends not only on the local concentration of ab- 
sorbers, but also on its filling factor (/), and, in lesser 
degree, both on the size of inhomogeneities, and on 
their topological properties. In the limit of small scale 
inhomogeneities, when the mean free path length in 
the far wing exceeds a typical size of the inhomogene- 
ity, the optical depth in the wing may be reduced to 
a simple form At = (ari/a;^)/|Aa;|, where ri is the 
Sobolev optical depth and /S.x — {dx/dl)Al is the fre- 
quency shift in comoving frame at the length interval 
Al. Under this approximation the Ha profile from 
the inhomogeneous core is the same as that from the 



homogeneous core with the "effective optical depth" 
Ti.cff = fTi, The analysis of the Ha damping wings 
thus permits us to estimate the Ha effective optical 
depth in the mixed core, or, equivalently, fn2 {112 be- 
ing the population of the second hydrogen level). It 
is noteworthy that the derived value of fn2 is inde- 
pendent of the adopted distance of a supernova. 

3. Identification of Ha damping virings effects 
in SN 1997D 

The spectra of the type II-P supernova SN 1997D 
with a low expansion velocity have been obtained at 
the early {t « 150 d) and the l ate {t « 300 d) neb- 
ular epoch by Turatto et al. (1998). The nebular 
spectrum on day 300 was reproduced in the model 
of a spherically symmetric mixed core with the to- 
tal ejecta mass of 6 Mq and the kinetic energy of 
10^" erg (Chugai & Utrobin 1999). However, it was 
stressed there that Ha on day 150 shows a significant 
redshift in a disparity with a spherically symmetric 
model. Note, that the asphericity of the ^^Ni distri- 
bution might produce the Ha asymmetry on day 150 
likewise it was in SN 1987A. However, this reason for 
SN 1997D is discarded by the absence of Ha asym- 
metry on day 300. We will see below that damping 
wings effects may resolve the revealed controversy. 

The assumed envelope model is similar to that of 
the Section 2, i.e., it consists of a line-emitting homo- 
geneous core and a scattering halo. For models con- 
sidered here the core velocity is vi = 650 km s~^ and 
the outer velocity of the scattering halo is V2 = 2000 
km s~^; the latter is consistent with the steep density 
drop in the range v > 1500 km for the hydrody- 
namical model of SN 1997D (Chugai & Utrobin |l999| ). 
The Sobolev optical depth in the halo is assumed to 
decrease outward as T2{1 — v/v2) with T2 = 10. An 
additional factor implemented into the model is the 
Thomson scattering, although it is of minor impor- 
tance. The corresponding optical depth (tt) of the 
mixed core, estimated from the Ha luminosity (dis- 
tance D = 13.43 Mpc) using data by Turatto et al. 
( |1998| ), is 0.16 on day 150 and 0.08 on day 300. Bel- 
low, I adopt Tt = 0.1 for t = 150 d and tt = 0.05 
for t = 300 d to roughly take into account that the 
hydrogen filling factor in the mixed core is essentially 
less than unity (Chugai & Utrobin 1999| ). 

In the late time spectrum on day 300 the model 
without damping wings is not very much different 
to the observed profile (Fig. 2). It implies that 
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the damping wings at this epoch are not significant, 
although, the model with the damping wings and 
Ti,efF = 5 X 10^ fits the observed profile better. On 
day 150 the model without damping wings is strik- 
ingly dissimilar to the observed profile: the latter is 
broader and shows both significant redshift and red 
wing. The model with damping wings included and 
''"i.eff = 5 X 10^ turns out to be successful in the de- 
scription of the observed profile. Note, that with the 
core velocity set by the late time profile, the primary 
fitting parameter is the effective optical depth in the 
mixed core. 

The found effective optical depth value being com- 
bined with the Sobolev optical depth may be used to 
derive the hydrogen filling factor in the mixed core. 
To estimate the Sobolev opt ical d epth I adopt 0.002 
O of '^^Ni (Turatto et al 

V 

g cm 



M, 



in the mixed core 

13 



< 650 km s-i). The density p « 5 x 10 
is adopted on day 150 fo llowin g the model 
of SN 1997D (Chugai & Utrobin |l99S| ). The non- 
thermal excitation of hydrogen may be calculated as- 
suming that all the energy of gamma-rays from the 
radioactive decay of ^^Co is uniformly deposited in 
the mixed core, which is a sound approximation at 
the early nebular epoch when the mean free path for 
gamma-rays is small. The energy deposited in H-rich 
matter is spent (apart from Coulomb heating), even- 
tually, on the excitation of the second level of hydro- 
gen. With the depopulation rate controlled by the 
two-photon decay and the coUisional de-excitation I 
find 712 ~ 1-3 X 10** cm^^ on day 150, which leads to 
the Sobolev optical depth of w 1.7 x 10^. The lat- 
ter being combined with the effective optical depth 
(Ti,off = 5 X 10^) results in the hydrogen filling factor 
of / « 0.3. The latter value is quite realistic, although 
it is a factor of 1.5 larger than the estimate fr om th e 
nebular model on day 300 (Chugai & Utrobin 1999 ). 
The difference is not dramatic taking into account 
possible uncertainties of both models. The factor of 
ten lower value of the effective optical depth in Ha 
on day 300 compared to that on day 150 is consis- 
tent with the decrease of the Sobolev optical depth in 
Ha expected for the mechanism of the non-thermal 
excitation (e.g., Xu et al. 1992 ). 

Although describing an inhomogeneity of the ex- 
cited hydrogen in the mixed core in terms of the filling 
factor seems sensible, a possibility exists that between 
the mixed core and unmixed hydrogen envelope there 
is a transition layer of strongly excited hydrogen with 
the filling factor of unity. To estimate qualitatively 



effect of this component I introduced in the model 
on day 150 an additional spherical layer in the veloc- 
ity range 650-700 km s~^ with filling factor of unity, 
negligibly small net emission, and the Sobolev optical 
depth being equal to the effective optical depth of the 
mixed core. Such a choice of optical depth suggests 
factor of three lower hydrogen excitation in the tran- 
sition layer than in the mixed core. The new model 
is found to produce the same Ha profile as in the 
old model provided the effective optical depth of the 
mixed core is of Ti^eff = 4 x 10^, i.e., 20% lower com- 
pared to the old model. This implies that the value of 
the Ha effective optical depth found from the profile 
modelling is quite robust to reasonable edditing the 
geometry of the distribution of the excited hydrogen 
in the mixed core. 

Summing up, the Ha damping wings effects pro- 
vide a natural explanation for the dramatic discrep- 
ance between the model based on the Sobolev approx- 
imation and the observed profile in SN 1997D on day 
150. This conclusion is strengthen by the absence of 
a reasonable alternative. 

4. Summary and Discussion 

At the early nebular epoch of SN II-P the Ha emis- 
sion should exhibit observational signatures of damp- 
ing wings. Although, the damping wings were never 
considered before as relevant for supernova optical 
spectra, the claimed effects are robust and may be 
easily identified. We see damping wings effects in 
SN 1997D spectrum, where they are responsible for 
the redshift, large width and extended red wing of 
Ha on day 150. The detection of Ha damping wings 
effects in SN II-P provides a direct estimate of the 
Ha effective optical depth (i.e. the product of the 
Sobolev optical depth and the hydrogen filling fac- 
tor) in a distance independent way. We thus obtain 
a reliable test for nebular models of SN II-P. On the 
other hand, with a confident nebular model in hand, 
we may able to measure the hydrogen filling factor 
in the mixed core. The price for this knowledge is 
the abandoning the Sobolev approximation and the 
need for the full radiation transfer treatment of the 
Ha radiation. 

Needless to say that the simple model of the Ha 
profile used above for SN 1997D admits some im- 
provements. For instance, one may take the radial 
dependence of hydrogen excitation and filling factor 
in the mixed core into account. Another modifica- 
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tion, might be the inchision of the absorption of Ha 
quanta in the Paschen continuum. This may be done 
if the radiation transfer and the hydrogen excitation 
are calculated simultaneously. More detailed descrip- 
tion of the hydrogen inhomogeneities is also welcome, 
although the expected change of the value of the ef- 
fective optical depth would be within 20%. 

To recover the damping wings in other SN II-P it 
is necessary to use at least two nebular spectra: at 
the early epoch, (around day 150) and at the late 
epoch, but before the dust formation, (e.g., between 
300 d and 400 d). The early nebular epoch is most ap- 
propriate to extract the information about damping 
wings, while the late time nebular spectrum provides 
a template to set the core velocity and to secure re- 
sults from a confusing effects of asymmetric ^^Ni dis- 
tribution. Notably, it is ^^Ni asymmetry that prevent 
us from using a simple symmetric model to study Ha 
damping wings effects in SN 1987A. 

I am grateful to Bruno Leibundgut for hospital- 
ity at ESO and to Massimo Turatto for the generous 
permission to use the spectra of SN 1997D. This work 
was supported by RFFI grant 98-02-16404. 
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Table 1: Parameters of demonstration models 
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Vl 


n X 10"* 


V2 


T2 




(km s-i) 




(km s-i) 




Ml 


2000 


0.001 


2000 





M2 


2000 
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2000 





M3 


2000 
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2000 
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2000 


0.001 


4000 
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M5 


2000 
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10 
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Radial velocity (1000 km s 

Fig. 1. — Effects of Ha damping wings expected for 
SN II-P. The panel (a) displays the case without con- 
tinua and a scattering halo. The dotted line is the 
Sobolcv approximation (model Ml), while the solid 
line (model M2) and the dashed line (model M3) show 
effects of damping wings. The panel (6) displays the 
case with the continua and the scattering halo. The 
dotted hne is the Sobolev approximation (model M4), 
while the solid line shows the model M5 with damping 
wings. 
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Fig. 2.— Calculated and observed Ha in SN 1997D 

on days 150 and 300. The dotted line is the observed 
spectrum; the Sobolev approximation is shown by the 
thin solid line, while the thick solid line shows the 
model with damping wings. 
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